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One of the key ingredients of the Unified Model of Active Galactic Nuclei (AGN) is the presence 
of a torus-like optically thick medium composed by dust and gas around the putative supermassive 
black hole. However, the structure, size and composition of this circumnuclear medium are still 
matter of debate. To this end, the search for column density variations through X-ray monitoring 
on different timescales (months, weeks and few days) is fundamental to constrain size, kinematics 
and location of the X-ray absorber(s). 

Here we describe our project of mining the Swift-X RT archive to assemble a sample of AGN 
with extreme column density variability and determining the physical properties of the X-ray 
absorber(s). We also present the results obtained from a daily-weekly Swift-X RT follow-up mon¬ 
itoring recently performed on one of the most interesting new candidates for variability discovered 
so far, Mrk 915. 
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1. Variable Swift- AGN sample 

According to the Unified Model of Active Galactic Nuclei (AGN; [1]), an obscuring optically 
thick medium composed by dust and gas arranged in a torus-like geometry is present around the 
nuclear engine. However, the structure, size and composition of this circumnuclear medium are still 
matter of debate; important constrains have been recently provided by the study of the absorption 
variability, almost ubiquitous in bright absorbed AGN ([2]). X-ray absorbing column densities, 
Ah, have been observed to change over different timescales, from a few years down to weeks or 
days, for a handful of AGN (e.g., [3] and references therein), allowing us to investigate the X-ray 
properties of the absorbers down to sub-parsec scales. The emerging picture is that multiple neutral 
and ionized absorbers co-exist around the central black hole, located at different distances from it. 

We recently started a project aimed at finding more examples of these sources in order to 
study the physical properties of the X-ray absorbers on a larger statistical basis. We exploited 
the vaste amount of X-ray observations provided by the X-ray Telescope (XRT; [4]) onboard the 
satellite Swift. XRT long exposure observations are always split in segments of ~ 1 ks, with typical 
integration times of 5ks per day. This generates, for bright sources, a natural short-timescale 
monitoring, with spectra good enough to perform a basic spectral analysis. As an example, for an 
AGN with a flux of fb. 5 -iokeV ~ 5 x 10~ 12 ergs cnr 2 s -1 , in ~ 5ks we expect to collect ~ 500 net 
photons. This makes the Swift-XRT archive an important database to search for different timescale 
variability of bright sources. To assemble our sample, we considered: 

1. the serendipitous sources detected in the field of view of XRT observations performed as 
follow-up of 76 gamma-ray bursts. These observations correspond to ~ 49 square degrees; 

2. the XRT sources classified as AGN that are possible counterparts of objects detected by the 
Burst Alert Telescope (BAT) onboard Swift [5, 6]. This AGN sample can be considered 
unbiased in term of column density, at least for Ah < 10 24 cm 2 (a limit above which a 
source falls in the so-called Compton-thick regime). In the energy band covered by BAT 
(14— 195 keV), the intrinsic emission is indeed weakly affected by absorption levels below 
few 10 24 ciW 2 . 

For each source, we split the XRT data sets on a daily basis, combining all the observation segments 
within 1 day (for the AGN-BAT sample) or 2 days (for the serendipitous sources). We filtered out 
time bins shorter than 200sec. The light curves span a wide range of number of time bins (from 
one up to more than 200 spectra) and timescales (from few days to years). For 299 serendipitous 
sources and 405 BAT AGN we have at least two spectra. 

As a proxy of the amount of column density we used the hardness ratio HR = (H — S) / {H + S ), 
where S = 2 — 4keV and H = 4— lOkeV. In the 2 — lOkeV energy band, the spectrum of an AGN 
can be approximate, at the first order, by an absorbed power law, while below ~ 2keV different 
soft components (soft excess, reflection, scattering, etc) are also present. Under this hypothesis, the 
HR value as defined here provides a strong indication of the amount of absorption: in Fig. 1 (left 
panel) we compare the HR vs. Ah derived from the analysis of the XRT data with the HR measured 
for a number of simulated spectra (photon index I' = 2) with different values of absorbing column 
density. For each input Ah value, we indicate with different colours the 16 th and 84 th percentile of 
the distribution of 100 tests with three different values of 2 — lOkeV source photons (from lightest 
to darkest, 20, 100 and 500). From the plot it is evident that the HR is a proxy of the Ah, and 
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/V H (onr 2 ) HR intrinsic scatter 

Figure 1: Left panel: HR of simulated spectra (T = 2) for different values of /V H . Shaded bands: 16 th and 84 th percentile 
from 100 simulations assuming source counts from 20 (yellow) to 500 (purple). Superimposed, we plot the real values 
derived by the XRT data for the BAT AGN considered here (black dots and error bars). Right panel: Examples of 
posterior probability of HR intrinsic scatter for 3 sources for which spectral variability has been observed at 99.99...% 
confidence level. 
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Figure 2: Observed HR light curves derived from the XRT archival data binned on a daily basis for two of the variable 
candidates (Mrk 915, left ; and ESO 323-077, right). 

variation in the HR can be interpreted as an indication of a change in the absorbing column density, 
at least for /V H between ~ 10 22 cm~ 2 and few times 10 23 cnr 2 . 

In order to probe the spectral variability, for each source we measured the possible non-zero 
amplitude of an intrinsic HR scatter superimposed to the statistical one. We adopted a Bayesian ap¬ 
proach similar to the one used by [7] to determine the intrinsic scatter of cluster velocity dispersion 
in presence of measurement errors, adapted here to Poisson-based quantities. Basically, starting 
from the source and background photons in the two bands (2 — 4keV and 4 — lOkeV) and in the 
different time bins, we calculate the probability for each value of intrinsic scatter given the data. 
As an example, the probability distributions for 3 sources for which a spectral variation is clearly 
detected are shown in Fig. 1 (right panel), where we plot the probability associated to an intrinsic 
scatter value ranging from 0 (corresponding to the no-scatter case) to 0.6. As evident from the 
figure, the probability of null intrinsic scatter is <C 10 3 , implying a negligible probability that the 
measured HR variability (shown in Fig. 2 for two of the sources) is due to statistical fluctuations. 

By applying a very conservative threshold of 10 3 in the probability of null intrinsic scatter, 
we found 15 AGN candidates for variable absorbers, 2 of which are serendipitous and 13 are 
BAT-selected. The sources are listed in Table 1, together with the main characteristics of the 
XRT archival data considered here. For the AGN observed by BAT, we summarise also the basic 
information extracted from the Swift- BAT 70 month survey catalogue 1 . Most of the BAT-selected 
sources are well-known variable AGN: this confirms the validity of the adopted method. We are 

'See http://swift.gsfc.nasa.gov/results/bs70mon/ 


3 





































AGN absorber structure through the XRT archive 


Lucia Ballo 


Source name 

XRT archive 

BAT catalogue 

Ref. 


Tot. Exp. 

(count rate) 

# spectra 

^!4-l95keV 

M4-195 keV 

z 

Opt. class. 


d) 

(2) 

(3) 

(4) 

(5) 

(6) 

(7) 

(8) 

(9) 

Mrk 915 

23.2 

0.14 

6 

31±6 

1.7T0.3 

0.02411 

Syl 


NGC 1365 

154.8 

0.14 

34 

64±4 

2.0T0.1 

0.00546 

Syl.8 

[8, 9, 10, 11, 12] 

ESO 362-18 

96.9 

0.53 

31 

49 ±4 

1.9±0.1 

0.01245 

Syl.5 

[13] 

CGCG 031-072 

10.7 

0.03 

3 

21±5 

1.9T0.3 

0.03305 

Syl 


NGC 4395 

243.6 

0.11 

190 

25+3 

2.0T0.2 

0.00106 

Syl.9 

[14, 12] 

ESO 323-077 

7.5 

0.09 

3 

33 

2.1 

0.01501 

Syl.2 

[15] 

Mrk 335 

299.2 

0.34 

182 

18±4 

2.3T0.3 

0.02578 

Syl.2 

[16, 17, 18, 3] 

NGC 3227 

14.5 

0.32 

6 

10±4 

2.0T0.1 

0.00386 

Syl.5 

[14, 3, 12] 

NGC 3516 

36.8 

0.46 

31 

118±4 

1.9T0.1 

0.00884 

Syl.5 

[19, 20, 3, 12] 

NGC 4051 

39.6 

0.67 

32 

40±3 

2.3T0.1 

0.00233 

Syl.5 

[21,22, 14,3, 12] 

Mrk 766 

108.5 

0.58 

56 

22±3 

2.5±0.2± 

0.01293 

Syl.5 

[23, 24, 25, 26, 14] 

NGC 3783 

21.7 

2.2 

6 

181±5 

2.0T0.1 

0.00973 

Syl 

[27, 14, 3, 12] 

MCG+11-11-032 

14.6 

0.02 

3 

18±4 

1.8T0.3 

0.036 

Sy2 


SWJ012118-125727 

53.6 

0.01 
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SWJ105826-231414 

59.4 

0.004 
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Table 1: Swift-based information for the 15 objects that our method candidates as variable with high confidence. We 
list in the first part the 13 BAT-selected AGN, and in the second part the 2 serendipitous sources. Col. (1): Source name. 
Cols. (2) - (4): Information related with the XRT archival data used in this work: total exposure time, in units of ksec; 
average 0.5 — lOkeV count rate; number of useful spectra. Cols. (5) - (8): Information derived from the Swift-BXT 
70 month survey catalogue: BAT flux in the 14— 195keV energy band, in units of 10~ 12 ergs cm - " s _1 ; spectral index 
(computed from a power-law fit to the 8-band BAT data); redshifts (taken from the online databases NED and SIMBAD); 
optical classification. Col. (9): Not-exhaustive list of previous studies of variability in the X-ray band. 

exploring its potentiality by relaxing the threshold in probability of null intrinsic scatter. However, 
even the more restrictive threshold used here demonstrates the effectiveness of this method: it 
provided us with 5 new candidates (two of which are serendipitous) that were not previously known 
as variable sources. Below we present the results obtained from a daily-weekly Swift-XRT follow¬ 
up monitoring that was recently performed on one of them, Mrk 915. 

2. A new candidate for variable absorber: Mrk 915 

Mrk 915 (z = 0.024) is a Seyfert 1.5-1.9 galaxy hosting a black hole of mass M bh = 0.6 — 
1.8 x 10 8 M q ([28]). Spectroscopic observations performed in the optical band ([28, 29]) suggest 
that part of the optical absoiption affecting the central regions originates outside the obscuring 
torus, and it is most likely associated with the dust lanes seen crossing the central source ([30]). 
For this source the HR variation observed has a probability less than 10 6 of having a null intrinsic 
scatter superimposed to the statistical fluctuations (Fig. 2, left panel). 

We were awarded 12 observations with XRT (covering ~ 3 weeks, for a total of 76 ksec net 
exposure time), to further investigate the Nu variability and its timescales. In Fig. 3 we show the 
temporal behaviour of Mrk 915, as reconstructed considering all the Swift datasets (archival on 
the left panel, new monitoring on the right panel). The data are binned both per observation (red 
filled circles) and per snapshot (grey dots). Comparing the 0.3 — lOkeV light curves from the 
new monitoring and the archival data, we confirm, on shorter timescales, the count rate variability 
previously detected (see Fig. 3, upper panels). However, the analysis of the ratio of count rates 
observed in hard and soft bands (H/S) highlights a change in the behaviour of the source: in the 
archival observations Mrk 915 underwent to a dramatic spectral variation on a timescale of 1 month. 
This is clearly revealed by the decreasing of the 4 — lOkeV vs. 0.3 — 4keV ratio of a factor of 3 
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Figure 3: Light curves of the 0.3 — lOkeV count rate ( upper panels ) and of the count rate ratios in the energy ranges 
2— lOkeV vs. 0.3 —2keV (middlepanels) and 4— lOkeV vs. 0.3 — 4keV (lowerpanels). We compare the trend observed 
in the archival data of Mrk 915 (left panels) with the results from the September 2014 Swift monitoring (right panels). 
As time zero-point, the starting time of the first archival observation is assumed. Grey dots: one bin per snapshot; red 
filled circles: one bin per observation. To highlight the observed variation in the H/S, we mark with grey dashed lines 
the mean ratio observed in the first observation of the monitoring. From Severgnini et al. (in prep.). 

(see Fig. 3, lower panels). A slightly less intense variation is also observed in the 2 — lOkeV 
vs. 0.3 — 2keV ratio (about a factor of 2.4). During the new Swift monitoring the change in the 
spectral shape is less prominent, of a maximum factor of ~ 1.5 on timescales of ~ 3.5 — 4days. No 
difference in the maximum amplitude of variation is observed when comparing the light curve of 
the 4 — lOkeV vs. 0.3 — 4keV ratio with the light curve of the 2 — lOkeV vs. 0.3 — 2keV ratio. 

To summarise, the Swift monitoring allowed us to find Mrk 915 in a different state than previ¬ 
ously detected, and to follow the source during a spectral variation. The origin of this behaviour is 
still under investigation; a detailed spectral analysis will be presented in Severgnini et al. (in prep.). 
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